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Abstract: Surface Acoustic Wave (SAW) technology is low cost, rugged, lightweight, extremely 
low power and can be used to develop passive wireless sensors. For these reasons, NASA is 
investigating the use of SAW technology for Integrated Vehicle Health Monitoring (IVHM) of 
aerospace structures. To facilitate rapid prototyping of passive SAW sensors for aerospace 
applications, SAW models have been developed. 

This paper reports on the comparison of three methods of modeling SAWs. The three models are 
the Impulse Response Method (a first order model), and two second order matrix methods; the 
conventional matrix approach, and a modified matrix approach that is extended to include internal 
finger reflections. The second order models are based upon matrices that were originally developed 
for analyzing microwave circuits using transmission line theory. Results from the models are 
presented with measured data from devices. 
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1. Introduction 

First order models of SAW devices are based upon the Impulse Response[l, 2], These models do 
not take into account second order effects such as internal reflections, frequency shifts, or allow for 
any physical arrangement other than equal electrode widths and spaces. For more accurate results, a 
matrix based approach was developed [3]. This approach has been further refined and modified to 
include internal finger reflections [4]. The reflections occur when the thickness of the metallization is 
sizeable enough to result in significant reflections. The extensions are based upon matrices that were 
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originally developed for analyzing microwave circuits using transmission line theory. The 
modifications are accomplished by breaking up the SAW device into zones, where the area under a 
metalized region is treated as one zone, and the area without metallization is treated as another zone. 
The impedance discontinuities that occur at the edges of the metal fingers enable the simulation of the 
internal reflections of the mechanical acoustic wave. The modifications also enable incorporation of 
the different velocities for each region, which produces a more accurate characterization of the 
frequency response of the device. The modifications also allow for unequal finger widths and spacing. 

In this paper, the three modeling methods are briefly discussed, then a comparison of the three 
models with results from prototype devices are presented. 


2. First Order Modeling using the Impulse Response Method 


The Impulse Response method [1] was used as the baseline for modeling the SAW device. This 
method is valid only for transducers where at least one of the two Inter-Digitated Transducers (IDTs) 
has a constant aperture or finger overlap [3]. This first order model includes both the mechanical and 
electrical behavior of SAW devices. It calculates the frequency response, the loss of the system, the 
admittance, and parameters for circuit simulators. This model assumes constant and equal spacing and 
finger widths. A simple circuit model (Fig. 1) can be used to convey the basic elements of the Impulse 
Response Model. The figure shows the source voltage and both the source and load impedances which 
are not part of the model. In the circuit model Ct is the total capacitance, B a (f) is the acoustic 
susceptance and G a (f) is the radiation conductance. 
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Fig. 1. (a) Basic SAW delay line and (b) the circuit model used in the Impulse Response Modeling. C T is the 
total capacitance, B a (f) is the acoustic susceptance, and G a (f) is the radiation conductance. 

The frequency response of a SAW device can be calculated by using the Impulse Response Model 
and is given by [5]: 


H{f) = 201og 




( 1 ) 


Where k is the piezoelectric coupling coefficient, C s is the capacitance per finger pair and unit 
length, H a is the aperture or overlap height of the fingers, fn is the center or synchronous frequency, N p 
is the number of finger pairs, /is the frequency, D is the delay length between the IDTs, and Z, is the 
wave impendence and X is given by: 


X = N p n 


(/-/o) 
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3. Conventional Matrix Method 


The matrices discussed in this work were originally developed for analyzing microwave circuits 
using transmission line theory. The methodology utilizing transmission matrices was modeled on the 
approach given by Campbell [3]. This method is based upon the Mason equivalent circuit using the 
crossed field technique (Fig. lb) [4], In this method, for modeling purposes, a SAW device is modeled 
as two IDTs with two electrical ports, and two acoustic ports. The acoustic ports represent mechanical 
waves traveling into and out of the IDT. The electrical ports represent the current and voltage of the 
IDT. The matrix for a SAW delay line is simply the multiplication of the matrices for the two IDTs 
and a matrix for the delay in between (Fig. 2). 
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Fig. 2. Transmission line matrix model of a complete SAW delay line. 


The transmission matrix is often used because it can be cascaded easily. The complete SAW device 
matrix is given by: 

[SAW{f)]=[T x (f)D,(f)T 2 (f)\ (3) 

Where / is the frequency, Ti(f) is the transmission matrix for an IDT, and D(f ) is the delay matrix. 
The delay matrix (D(f)) is modeled after an acoustic transmission line and is given by: 

2f* d 1 

V 0 

’ (4) 

0 e 1 v 

where X is the wavelength at the center or synchronous frequency and d is the distance between the 
reference planes, or in this case the center of the two IDTs. 
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Each IDT is modeled separately as a single transmission matrix (Tj(f)) with one electrical port, and 
two acoustic ports (Fig. 3). 
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Fig. 3. Transmission matrix model of an IDT. 

This allows the acoustic waves (W,) and electrical parameters (a; and b;) to be related through the 
use of transmission matrix Ti(f) in: 


( 5 ) 
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The transmission matrix is decomposed into sub-elements, given by 
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The elements of the transmission matrix are given by: 
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The symmetry factor s is based upon the number of fingers (N t ): 



5 = -\ N ‘ . 

(14) 

The transit angles 0 e , 0 C , 

and 9,, are given by: 



9 e (f) = (coC t +B a )(R s +Z e ) , 

(15) 


0 c (f) = oC,(R s +Z e ) , 

(16) 


O t (f) = N±S. 

(17) 

The combined load or source impedance is Z e , and R s is the combined IDT metal 
resistance. The frequency detuning parameter (<5) is given by 

and lead 


s Jn(f~Q_ kn _ 

V 

(18) 

where v is the velocity, and kn is the mutual-coupling coefficient and is given by: 



h ( h V 

Ki=Ki P + K\m^ + Ku t • 

A \A y 

(19) 


where h is the height of the metal fingers. The mutual coupling coefficient parameters k IIp , k IIm , 
and kus, are material dependent, for quartz they are 0.0004, 0.02, and 7.9 [3], The total capacitance for 
the IDT is given by: 


office.. 


( 20 ) 


where C s is the capacitance per unit length for a pair of fingers. The conductance ( G a ) is given by [3]: 


G a = %k 2 C fJN -If 


sin 






e. 


( 21 ) 


And finally, the susceptance ( B a ) is given by [3]: 

B a =\6k 2 C s m-\) 2 


v ^ y 


( 22 ) 
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4. Modified Matrix Method 


For more accurate results than are given by the first order approach or the conventional method, the 
matrix approach was extended to include second order effects such as the internal finger reflections [3, 
6]. The model divides an IDT into !4 wavelength sections. These sections are further divided into 
zones. Two of the zones are un-metalized areas (1/8 of a wavelength) around one zone that is 
comprised of a metal finger (1/4 of a wavelength). Each zone is modeled by a transmission line matrix 
equivalent circuit (Fig. 4). Two identical circuits model the un-metalized areas, while the middle 
circuit models the area under the metal finger. The transmission matrix relates the SAW voltages V i 
and V 2 , to the currents 1 1 to I 2 respectively. The acoustic wave is assumed to have entered from the left 
and travels through the element towards the right. In this model Z u and Z m are the acoustic 
impedances for the metalized and un-metalized areas, C 0 is the capacitance for a single finger, 0 U and 
0 m are the acoustic angles of the substrate, and the turns ratio of the transformer is assumed to be 1:1. 


A A A 



Fig. 4. Crossfield Model using Mason Equivalent Circuit for a V 2 wavelength section of an IDT. 

The transmission matrix representing the middle area shown in Fig. 4 is the circuit for a metalized 
region that is assumed to be lossless and is given by: 

coshO 0 m (f)) Z m sinh(//? m (/)) 

-Uinhij^QO) cosh(//U/)) ' (23) 

Am 

The transmission matrix (23) is determined by the acoustic angle 0 m and the metalized region’s 
impedance Z m . The impedance Z m is calculated with: 



z„Xf) 


k 2 C s H a f m ’ 


(24) 
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where k 2 is the piezoelectric coefficient, C s is the capacitance for a single pair of electrodes per unit 
length, and H a is the aperture width. The acoustic angle of the substrate 0 m , is given by 

W) = 7 7 - > (25) 

where / is the frequency, and f m is the frequency shift caused by the velocity change as the acoustic 
wave travels under the metalized area: 


L = V -J • (26) 

Where v m is the acoustic wave velocity under the metalized area and X is the wavelength of the 
device. The metalized velocity (v,„) is 3134 m/s for ST cut Quartz. 

The matrix (23) calculates the parameters for the metalized area, but cannot be used for the un- 
metalized sections. This leads to the transmission matrix ( R u (j )) for the un-metalized region as is given 
by: 




cosh(/T/(/)) 

jsttuvm 


2 , stah (jejf)) 

cosh (jO„(f)) 


(27) 


The transmission matrix (27) is determined by the acoustic angle 
impedance Z„. The impedance Z„ is calculated with: 
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z u (f) = 


1 

k 2 C,H a f 0 ’ 


(28) 


where fo is the synchronous frequency of the acoustic wave under for the un-metalized area. The 
acoustic angle of the substrate 0 U , is given by 


W) = 


x_f_ 

4/o ' 


(29) 


where fo is the synchronous frequency of the acoustic wave, which is calculated using the acoustic 
wave velocity under the un-metalized area. The un-metalized velocity (v) is 3158 m/s for ST cut 
Quartz. 

To find the transmission matrix for the 14 wavelength periodic element ( R T ( f)) one must multiply the 
three matrices together for both metalized region and the un-metalized regions adjacent to it: 

K(/)]=K(/)]K(/)]K(/)] • (30) 

To find the transmission matrix ( Q(F )) for an entire IDT one simply raises the (Riff)) matrix to the 
power of 2N P \ 


IDT X (/ )] = \R,(f)^ Np . 


(31) 


Where N p is the number of electrode pairs, so 2N P is the total number of electrodes in the IDT. 


7 



The matrix for a SAW delay line is simply the multiplication of the matrices for the two IDTs and 
the delay or space between the IDTs. The SAW matrix is given by: 


[iMW'O')] = [ZDr, (y)][A (/)][707- 2 (/)] . (32) 


5. Prototype Design 

A simple SAW delay line that consists of two identical un-apodized IDTs was chosen as a 
prototype to illustrate the capability of the models. Each IDT has 63 finger pairs that are 17 pm wide. 
The spacing between the fingers is 17 pm also. The center or synchronous frequency is 46.44 MHz, or 
a wavelength of 68 pm. The aperture height is 2730 pm. The delay length between the IDTs is 10 
wavelengths or 680 pm. The design was fabricated on two different ST cut quartz wafers. One with a 
single side polished, and one wafer with both sides polished. The aluminum thickness is 50 nm for the 
wafer with a single side polished and 250 nm for the wafer with both sides polished. Also, a 1 12.78 
MHz device was fabricated on the single side polish wafer. It has 90 pairs of 7 pm fingers, a wavelength 
of 28 pm, and an aperture height of 1 105 pm. The delay length is 104 wavelengths or 2.904 mm. 

6. Results 

All three methods adequately model the frequency response and amplitude for the main lobe and 
the first two side lobes for cases without any mass loading due to the metal fingers (Fig. 5). The 
modified matrix more accurately captures the frequency shift due to the mass loading of the metal 
fingers. In Figure 6, the ideal first order model (impulse response model) and the conventional matrix 
results are both centered about the synchronous frequency. The measured results and the Modified 
matrix results are both shifted down in frequency due to velocity changes from mass loading effects. 

When comparing the measured data from the two figures, it is noticeable that the main lobe peak of 
Fig. 5 does not have the same artifacts as are seen on the peak of the main lobe in Fig. 6. These 
artifacts are due in part from bulk waves that are reflected from the polished bottom surface of the 
wafer. The roughness of the non-polished surface disperses the bulk waves which results in 
diminished artifacts in the main lobe peak (Fig. 5). Also note that the peak in Fig 6 is shifted farther 
away by mass loading than the peak in Fig. 5. 



Frequency (MHz) 

Fig. 5. Comparison of model results with data from a single side polish wafer, with 50 nm of aluminum. 



Frequency (Mfiz) 

Fig. 6. Comparison of model results with data from a double side polish wafer, with 250 mu of aluminum. 

While the models work well for the side lobes of low frequency devices at 46.4 MHz, as the frequency 
increases the accuracy goes down for the side lobes (Fig. 7). The measured data shows higher order 
effects and noise at 112.8 MHz that is not present in the lower frequency device data. These results 
show that opportunities exist for improvements to both the models and the test setup. The inclusion of 
higher order terms to the models will increase their fidelity. The fidelity of the test setup can be 
improved by dicing the wafer or by the inclusion of gating functions in the time domain. Time gating 
will allow the removal of unwanted reflections from adjacent devices on the test wafer without having 
to dice the wafer and add absorber material to the devices. Essentially the addition of the gating 
functions will allow us to take data on the wafer, while achieving the results as if the devices had been 
packaged in an ideal manner. 



Frequency (MHz) 

Fig. 7. Comparison of model results with data from a single side polish wafer, with 50 nm of aluminum. 
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7. Conclusion 


In this paper, results from prototype devices are compared to the Impulse Response model, a 
conventional matrix model, and a modified matrix model. In the absence of mass loading from the 
metal fingers and the associated frequency shift, all three methods model the main lobe and next two 
side lobes fairly well. However, the results show that the modified matrix methods more accurately 
modeled second order effects such as frequency shift due to the metal thickness of the IDT fingers. 
Neither the Impulse response method nor the conventional matrix method model the frequency shift. 
Therefore, the modified matrix method is the most accurate for a wider range of parameters such as 
finger metal thickness. Future work will include a comparison the Coupling of Modes, the Impulse 
Response model, and the modified matrix model. Also, to achieve higher fidelity our future models 
will include higher order terms. To improve our wafer level test setup, time gating in the time domain 
will be added to remove artifacts that would not be present if the devices were packaged in an ideal 
manner with absorber materials. 
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